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Abstract 

By using the Rosenfeld density functional we determine the number density 

profiles of both components of binary hard-sphere fluids close to a planar hard 

wall as well as the corresponding excess coverage and surface tension. The 

comparison with published simulation data demonstrates that the Rosenfeld 

functional, both its original version and sophistications thereof, is superior to 

previous approaches and exhibits the same excellent accuracy as known from 

studies of the corresponding one-component system. 
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I. INTRODUCTION 



A variety of experimental techniques has emerged which allow one to resolve the inhomo- 
geneous density distributions of fluids at interfaces, a subject which enjoys broad scientific 
interest. In this context the ability to manufacture highly monodisperse colloidal suspen- 
sions has turned out to be particular useful as they provide the possibility to tune the 
effective interactions in these systems such that, e.g., the colloidal particles closely resem- 
ble hard-sphere fluids [[!]] . Since many of these experimental probes are indirect scattering 
techniques there is a substantial demand to guide them theoretically. Computer simulations 
and integral theories are important tools of statistical physics to address these issues. 
Density functional theory (DFT) || has emerged as an additional approach which is capa- 
ble to capture interfacial phase transitions and to sweep the thermodynamic and interaction 
parameter space of the system under consideration. The potential to combine these two 
possibilities poses already a major challenge for the other techniques. If DFT acquires in 
addition the same accuracy as the other two techniques, it could gain a clear competitive 
edge. 

Although there is no recipe for constructing systematically a reliable DFT in spatial di- 
mensions d > 2, the constant flux of developments over many years has led to a rather high 
level of sophistication. Among these theories for hard-sphere fluids, which act as paradig- 
matic systems and stepping-stones for more complicated models, the Rosenfeld functional 
has emerged as a particular powerful theory which resorts to the fundamental geometrical 
measures of the individual sphere Q. For the standard test case of the highly inhomoge- 
neous density distribution of a one-component hard-sphere fluid near a planar hard wall, 
the predictions of the Rosenfeld functional are very close to those of numerical simulations 
serving as benchmarks. For this case the mean square deviations [see, c.f., Eq. (p2|)j of 
Rosenfeld DFT results from the simulation data from Ref. |J are at most 1 x 10~ 3 at high 
packing fractions, otherwise less than 3 x 10 -4 . 

Another virtue of the Rosenfeld functional is that is easily lends itself to the generalization 
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to multi- component hard-sphere fluids. This opens the door to investigate rich new physical 
phenomena as particles of different size compete for interfacial positions ||. Even for the 
simplest multi-component system, the binary hard-sphere fluid, there are only relatively few 
theoretical studies which determine their structural properties near a planar hard wall, using 
Monte Carlo simulations J7|— PT0|] , integral equation theories [[J, and various kinds of density 



functional theory [§-[l5|], as well as in spherical pores [|l(|. Here we analyze this problem 
by using the corresponding Rosenfeld functional both in its original version f3| as well as 
for sophistications thereof [|17],[]1|]. By comparing these results with published simulation 
data [Kj we assess to which extent the quantitative reliability of the Rosenfeld functional 



for the one-component hard-sphere fluid remains valid for the corresponding binary system. 
Moreover we determine concentration profiles, the excess coverage, and the surface tension 
of the binary hard-sphere fluids at a hard wall. We describe the DFT in Sec. fl] and report 
our results in Sec. |T| followed by a summary and our conclusions in Sec. |V[ The Appendix 
contains important technical details. 

II. DENSITY FUNCTIONAL THEORY 

The Rosenfeld functional for the excess (over the ideal gas) Helmholtz free energy of a 
mixture of hard spheres with number density profiles {pj(r)}, i = 1, . . . , N, can be written 
as i 

= / d 3 r $({n Q (r)}) (1) 
which is a functional of the four scalar weighted densities n a (r) for the A-component mixture 



N 

n a (r)=J2 dV^rVf'tr-r'), a = 0, . . . , 3, (2) 
i=i ^ 



with AN scalar weight functions and two three-component vector weighted densities 



n Q (rj 

N 



n «« = E / dVp^rOu^r-r'), a = 1,2, (3) 



with 2N vector weight functions u}\ . The weight functions contain only information about 
the fundamental geometrical measures of a single sphere of species i, namely its volume, 
surface area, and radius Ri, i.e., in particular they are independent of the density profiles. 
The explicit expressions for the weight functions are given in the Appendix. Q({n a }) = 
$1 + $ 2 + $3 is a function of the weighted densities with M 



$x = -n log(l-n 3 ), (4) 

'-2 ~ ni ■ 
1 -n 3 



n-ino — ni ■ n 2 

$ 2 = JlJ — ?_a ( 5 ) 



and 



- n 2 n 2 • n 2 _ n 2 , 2 , 



*■= Ml-n,)» ^ 24,(1 i^ "-^ < 6) 

where £(r) = |n 2 (r)|/n 2 (r), which is the ratio of the modulus of the vector weighted density 
n 2 (r) and the scalar weighted density n 2 (r). We want to note that £(r) = in the bulk and 
is small for small inhomogeneities. While this original Rosenfeld functional describes very 



successfully the fluid phase of a one-component hard-sphere system [Bl, it fails to predict 



the freezing transition. This failure has been studied in detail in Refs. []I7| and ||18|| . For the 
freezing transition it turns out that the zero-dimensional limit of the functional, in which a 
small cavity can accommodate only a single sphere, plays a key role. In a crystal the thermal 
vibrations around a lattice site can be interpreted as the motions in such a cavity formed 
by the neighboring spheres. Only if the statistical mechanics in such a cavity is described 
properly by the density functional, the freezing transition is predicted correctly. This is 
not the case for the original Rosenfeld functional. This problem can be fixed by modifying 
slightly the contribution $ 3 [see Eq. (§)] such that the freezing transition is predicted by the 
modified functional while at the same time, for lower packing fractions, the accuracy of the 
original functional in describing the inhomogeneous fluid is kept. The following modifications 



have been suggested [[T^HI 



*3. = ^(l-« 2 )' (?) 



with q > 2 and 



^ nt = ^(i-se+2e). (8) 



The first suggestion, $3, g , is an antisymmetrized version of $3 in Eq. (^) and the second, 
$3,1,14, interpolates between $ 3 of Eq. (|6|) and $3,00 i n the exact zero- dimensional limit 

While the modified Rosenfeld function with $ 3j od does successfully predict the freezing 
transition of the one-component system, it leads to modified bulk properties and hence 
cannot describe the hard-sphere fluid as accurate as the original Rosenfeld functional. We 
note that the difference between $3 of the original Rosenfeld functional and both $ 3 (J with 
q = 3 and $3,^ is of the order of 0(£ 3 ). Therefore we expect the biggest differences between 
the various versions of the Rosenfeld DFT to occur close to the wall where £ is largest. 

Both the original Rosenfeld functional and the modifications corresponding to Eqs. (|7|) 
and (|8|), i.e., the functionals that share common bulk properties, are very successful and 
accurate for the one-component fluid. But far less is known for binary mixtures. While in 



this latter respect in Ref. [20|] very good agreement between the density profiles obtained 



from the Rosenfeld functional and simulation data from Ref. f?j] has been mentioned, in a 



recent study JTO ] significant deviations between the Rosenfeld DFT results and simulations 



have been found. In this latter study an alternative but equivalent formulation pl| , |22[| of 
the original Rosenfeld functional has been applied. 

Here we are interested in the equilibrium density profiles p s ,o( r ) an d P&,o( r ) of both the 
small and 6ig components of binary hard-sphere mixtures close to a planar hard wall. To 
this end we freely minimize the functional 

n[p s (r),p b (r)]=T[p s (r),p b (r)} + J2 [ dV Pi {v') (VUr') - (10) 

i=s,b^ 

which is written in terms of the functional 

F\Pa(*),Pb{*)] = Fid[p s (*), Pb(r)] + J^ ea: [p s (r), p b (r)] (11) 



with the exactly known ideal gas contribution jF id , 

m d = £/ dv Pi {v>) (H\i Pi {v>)) - 1) , (12) 

i=s,b 

with Aj the thermal wave length of species i. For the equilibrium density profiles pifi(r), 
i = s,b, the functionals T and Q reduce to the Helmholtz free energy and the grand canonical 
potential of the mixture, respectively; fi s and fib are the chemical potentials of the two 
species. The external potentials entering into Eq. (10) model the planar hard wall at z = 0: 

{oo, z < Ri, 
(13) 
0, otherwise, 

% — s,b, with z the normal distance from the wall. The external potentials prevent the 
centers of spheres of species i to approach the wall, located at z — 0, closer than Ri in which 
case they are in contact. 

In the absence of spontaneous symmetry breaking due to freezing, which we do not 
consider here, the profiles pi } o(z), i — s,b, depend only on the normal coordinate z which 
simplifies the minimization of the functional. 

Far away from the wall, i.e., in the bulk system, the vector weighted densities rix and 
n 2 and thus £ vanish. In this limit both the original Rosenfeld functional and the two 
modifications corresponding to Eqs. fl7|) and (§) reduce to the same bulk expression given 
by 

n-lUo Tin 

$ M k = -no In 1 - n 3 ) + + — — * 14 

1 — n 3 247r(l — n z y 

and hence they share the same bulk properties. We want to emphasize that as a consequence 
of this feature all versions of the Rosenfeld functional predict density profiles which show 
the same asymptotic decay towards the bulk value ||23| . The weighted densities in the bulk 
limit are obtained by inserting the bulk densities p^bulk '■= Pi,o( z = oo) into Eq. @ yielding 

n 3 = R ^PiMik =^Vi, (15) 

i=s,b i=s,b 



n 2 = 4:71 ^ RfPiMik, (16) 

i=s,b 

TT-l = RiPi,bulki (17) 

i=s,b 

and 

n o=^2 PiMik- (18) 



The equation of state following from Eq. (|H|), 



Pp= ^- + ^^.L^i (19) 

l-n 3 (l-n 3 ) 2 127r(l-n 3 ) 3 



is the Percus-Yevick compressibility equation of state of the mixture [24|. This expression 



is related to the contact values of the density profiles according to the sum rule [25] 



pp= j2p^ z=R ^+°y ( 2 °) 



This sum rule is respected by the Rosenfeld functional as by any weighted-density DFT [26 
and therefore provides a test for the numerical accuracy of the calculations. In the following 
we suppress the subscript which indicates equilibrium profiles as opposed to variational 
functions. 

III. STRUCTURAL AND THERMODYNAMIC PROPERTIES 



A. Density profiles 

The number density profiles of both components of binary hard-sphere mixtures close to 
a planar hard wall are obtained by a free minimization of the functional given in Eq. (|10|). 
We use the original Rosenfeld functional as well as the modified versions corresponding to 
Eq. (§) and to Eq. (|7|) with q = 2 and q = 3. The systems considered here have two different 
size ratios, Rf, : R s = 5 : 3 and i?& : R s = 3 : 1, and various packing fractions t] s and r)b of the 
small and big spheres j2"7j| , respectively. The resulting density profiles are compared with 
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simulation data published in Ref. jn|. In addition we calculate the local concentrations 



& s ( z ) an -d ®b( z ) °f t ne small and big spheres, respectively, defined as 

*i(z)= (f+fv * = *,*>• (21) 
Ps {z) + p h {z) 

We find excellent agreement between the density profiles of both components obtained 
by density functional theory and the simulation data for all systems under consideration. 
This holds for all versions of the Rosenfeld functional analyzed here. While at low total 
packing fractions rj = r] s + rjf, the density functional theory results for all versions of the 
Rosenfeld functional are practically equivalent, small deviations among the results from 
different versions of the functional appear at larger values of rj, i.e., for t] > 0.3. We quantify 
the degree of agreement between our DFT density profiles Pi(z) and the simulation data from 



Ref. fLOl , available as data points (zj, pf m (zj)), j = 1 . . . Nf im and i = s,b, by determining 



the mean square deviations E iy i — s,b, defined as 

(22) 



1 ATsim / j \ 
1 .7 = 1 V 



' pf m (zj) - Pi ( Zj y - 



Pi,bulk 



.r- 

We find that E s and Ef, are at most 5 x 10~ 4 and 6 x 10~ 3 , respectively, for all versions 
of the Rosenfeld functional. However, since the statistical errors in the simulation data are 
comparable with or even larger than the differences between the density profiles obtained 
by different version of the Rosenfeld DFT this approach does not enable us to determine 
which of the various versions is the most accurate one. 

To illustrate the agreement between the DFT and the simulation data, in Fig. [I] we show 
the density profiles of the small spheres (a) and of the big spheres (b) for rj s = 0.0607 and 
776 = 0.3105 and a size ratio : R s = 5 : 3. The symbols (□) denote the simulation data 



from Ref. [10] and the solid line is obtained by the original Rosenfeld functional. The dotted 



lines denote coarse grained densities pf \ i = s,b and j = 0,1, defined as 



^(i) 



1 f j_ Pi( z ) 



z u+v_ z u) j PiMlk 



dz (23) 



.0) 



with zf^ = and the position of the jth minimum ||28|| . All details of the density 

profiles found in the simulations are reproduced very accurately by the density functional 
theory. The oscillatory behavior, i.e., the amplitudes, the phases, and the decay of the 
oscillations as obtained by DFT agree excellently with the simulations. The total packing 
fraction of the system r\ = T] s + rjb = 0.3712 is already rather high, giving rise to the 
pronounced structure of the density profiles. The corresponding concentration profiles & s (z) 
and <f>b{z) of the small and big spheres, respectively, are shown in Fig. These concentration 
profiles demonstrate that, apart from the purely geometrical constraints, near the wall the 
big particles are enriched and the small particles depleted. This is in line with the expectation 
based on the attractive depletion potential near a hard wall of a single big sphere immersed 
in a fluid of small spheres ||30|| . Correlation effects reverse this relative distribution in the 
second layer and restore it in the third. 

As mentioned above, small differences between the DFT results corresponding to the 
various versions of the Rosenfeld functional can be found for these values of 77. In order to 
be able to resolve these small differences magnified parts of the density profiles from Fig. [I] 
are shown in Fig. |3] together with the simulation data from Ref. ||10||(D). The solid lines in 
Fig. [3] correspond to the original Rosenfeld functional, the dotted lines correspond to the 
interpolated version [Eq. (|8|)] whereas the dashed and dashed-dotted lines correspond to the 
antisymmetrized version [Eq. (0)] with q = 2 and q = 3, respectively. All DFT results are 
very close to the simulation data. However, the deviations between the simulations and the 
antisymmetrized functional with q = 2 seem to be systematically the largest. 

In Figs. |4] and § we show the density profiles of a binary mixture with size ratio : 
R s — 3 : 1. The packing fraction of the small spheres is r] s = 0.0047 and that of the big 
spheres is 77^ = 0.3859 so that the total packing fraction rj = i] s +7]^ = 0.3906 is again rather 
high. Therefore there is a strong spatial variation of the density profiles. The agreement 
between DFT (solid line) and simulations (□) is again found to be excellent for both the 
density profile of the small spheres (a) and of the big spheres (b). The dotted lines denote 



9 



the coarse grained densities as defined in Eq. fl23|) ||29|| . The concentrations profiles of the 



small and big spheres, corresponding to these density profiles are shown in Fig. []. For this 
larger size ratio the anti-correlated behavior of & s {z) and $b (z) is even more pronounced 
than for the smaller ratio discussed in Fig. |2] and strongly locked in without additional fine 
structure such as the double peak appearing in Fig. 

In addition we test the numerical accuracy of our calculations by means of the sum rule 
given in Eq. (p0|) . In Table | the sum of the contact values of the binary mixture with 
size ratio Rb : R s = 5 : 3 for various packing fractions is compared with two equations 
of state. Pppy denotes the Percus-Yevick compressibility equation of state [Eq. ([19])], to 
which the Rosenfeld functional reduces in the bulk limit, and Ppmcsl corresponds to the 
more accurate Mansoori-Carnahan-Starling-Leland equation of state [31], which represents a 



generalization to a mixture of the very accurate Carnahan-Starling equation of state |32| for 
the one-component fluid. The very good agreement between the contact values and (3pp Y 
demonstrates the high accuracy of our numerical procedure. However, at higher packing 
fractions, j3p c PY deviates from the more accurate equation of state Ppmcsl- The same 
analysis of our results for a binary mixture with size ratio Rf, : R s = 3 : 1 is summarized in 
Table |TJ. 

Equation (^) represents a sum rule which must be fulfilled by the density profiles as 
obtained by any of the density functionals considered here. However, no corresponding rules 
are available for the individual contact values. We find that for all versions of the Rosenfeld 
functional under consideration here, the sum rule is respected equally well. However, the 
individual contact values may differ. This statement is in line with the expectation that 
the biggest differences between the various versions of the Rosenfeld functional occur in a 
region where £ is large, i.e., close to the wall, and is substantiated in Table |TJ for the binary 
mixture with size ratio Rb : R s = 5 . 3 and in Table [TV] for the size ratio Rb : R s = 3 : 1. 

Vested with this confidence in our numerical procedures we are now able to comment on 
the comparison of our DFT results for the original Rosenfeld functional with those obtained 



by earlier DFT studies |8Hl5| . We find that DFT approaches other than the Rosenfeld DFT 
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predict the structure of the density profiles of a binary hard-sphere mixture near a planar 
hard wall only qualitatively PJ9j,pJ]-p!5| , whereas the Rosenfeld functional yields quantita- 
tively reliable predictions. This is demonstrated by Figs. [I] and [3|-^ and the very small 
values of the mean square deviations E s and Ef, [Eq. (|22|)1. The predictions of the Rosenfeld 



functional agree excellently in all details with the simulation results from Ref. [IT]]. The 
deviations between the DFT results of Ref. |lOj, calculated with the Rosenfeld functional, 
and their own simulations originate most likely from numerically problems of the iteration 



procedure used in Ref. jlOfl . This suspicion is further supported by the fact that the DFT 



density profiles shown in Ref. |10| seem to violate the sum rule Eq. ([?(]). Thus we are led 
to the conclusion that the deviations between the Rosenfeld DFT results and the simula- 
tion data reported in Ref. |1(J most likely are artifacts generated by numerical problems 



in implementing the iteration procedure used in Ref. [|Hj] for solving the Euler-Lagrange 
equations corresponding to the Rosenfeld density functional. Therefore the doubts raised 
in Ref. [T(J about the performance of the Rosenfeld DFT for binary hard-sphere mixtures 
are not justified. We conclude that the Rosenfeld DFT exhibits the same high accuracy in 
predicting density distributions of binary hard-sphere mixtures as for the one-component 
hard- sphere fluid. 

B. Excess adsorption and surface tension 

One of the virtues of DFT is that based on the knowledge of the local structural properties 
Pi(z), i = s,b, it is straightforward to calculate also thermodynamic properties such as 
the excess adsorptions Tj and the surface tension 7. Here we determine these quantities 
near a hard wall for a binary hard-sphere fluid whose components exhibits a size ratio 
Rb : R s — 3 : 1. Our analysis is confined to the fluid phase of the mixture; the phase 
boundary for freezing is estimated from the bulk phase diagrams presented in Refs. |33| . 

The excess adsorption of species i, i — s, b, is defined as 
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dz (pi(z) - p iM ik) (24) 



o 

oo 

■s 

Ti/2 



This definition of the excess adsorption differs from the definition used in Ref. [33] as well 
as from that used in Ref. J35|. To recover the results for the excess adsorption of a one- 
component hard-sphere fluid in Ref. [34[] and Ref. |35j one has to subtract from and add to our 
results, respectively, the constant ^PiMik- These differences originate from different choices 
for the position of the wall. While for the one-component fluid there is no preference for 
any of these definitions, our choice used here appears to be particularly suited for a mixture 
because independent of the diameter Oi of species i the integrals in Eq. (|4j) start at the same 
lower bound for all species, namely at the position of the physical wall. We use the same 
definition for the position of the wall to determine the surface tension. Because the excess 
adsorptions follow from integrating over oscillatory functions, they depend very sensitively on 
the precise structure of the density profiles and require very accurate calculations. Moreover, 
near the phase boundary for freezing the original Rosenfeld functional yields values for T s 
and Tfe which differ significantly from those obtained from the modifications of the Rosenfeld 
functional. Thus, unless stated otherwise, we have determined the excess adsorption by using 
the modified Rosenfeld functional corresponding to Eq. (§), which is known to capture the 
freezing transition of the one-component hard-sphere fluid. 

In Fig. |7| we show the excess adsorptions T s of the small spheres as function of the packing 
fractions t] s and rjf,; note that r s (?7 s = 0,-%) = 0. For a fixed packing fraction of the small 
spheres rj s , the excess adsorption of the small spheres increases upon increasing rji,. The 
reason for this is, that the increasing packing effects of the big spheres associated with large 
values of % enforces also the packing of the small spheres, giving rise to a strongly enhanced 
contact value and very pronounced structures in the density profile of the small spheres. For 
very small rj s , i.e., r] s < 0.1, and large rjb we find that T s can become positive. As function 
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of rjb, T s exhibits a turning point for any fixed value of 77^ . 

The excess adsorption of the big spheres is shown in Fig. |^. For constant rj s and 
increasing rjb, Tb decreases. But due to the same mechanism as described above, r& increases 
upon increasing rj s for constant packing fractions of the big spheres rjb. The square symbols 
(□) in Fig. ^ denote simulation data for the excess adsorption of a one-component hard- 
sphere fluid near a hard wall at packing fractions 7^=0.3680, 0.4103, and 0.4364, respectively, 



taken from Ref. |35| . Whereas the simulation data for the two smaller packing fractions agree 
very well with our DFT results, that for the highest value of rj s clearly deviates from the DFT 
prediction. Thus from this comparison it remains unclear to which extent the simulation 
data and the DFT results are reliable close to freezing. 

In order to illustrate the large differences between the excess adsorptions calculated 
from different functionals, in Fig. || we show the excess adsorption T s and Tb calculated 
from the original Rosenfeld functional together with those calculated by using the modified 
functional corresponding to Eq. (Q). The packing fraction of the small spheres is rj s = 0.15. 
Only for small packing fractions of the big spheres, i.e., far away from freezing, we find 
good agreement between the results obtained from the different functionals. Close to the 
phase boundaries (BBl the differences become large, which is indeed surprising at first glance. 



However, the main contribution of the excess adsorption stems from the vicinity of the wall 
where £ is large and the different versions of the Rosenfeld functional are expected to differ. 
These differences were already indicated in the previous subsection by the behavior of the 
individual contact values afpi(Ri + 0), i = s, b, of the density profiles of the small and big 
spheres, respectively. For large distances from the wall the density profiles exhibit decaying 
oscillations and therefore are not expected to contribute essentially to the excess adsorption. 
Moreover, all versions of the Rosenfeld functional display a common characteristic decay 
because they share the same bulk properties so that the contributions to the excess 
adsorption far away from the wall are very similar for the original Rosenfeld functional and 
its modifications. 

The grand potential Q of a system in contact with a wall 
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^ — Qbulk + Qsurf (25) 

decomposes into a bulk contribution, Qbulk = —pV, given by the bulk pressure p in the system 
times the volume V occupied by fluid particles, and a surface contribution, Q sur f = jA, 
which is the surface tension 7 times the surface area A of the wall. Scaled particle theory 



(SPT) provides an approximate expression for 7 for a one-component hard-sphere fluid |J7] 



as well as a generalization to hard-sphere mixtures |24[ close to a planar hard wall. The 
surface tension of a one-component hard-sphere fluid within SPT is well tested and turns 
out to provide reliable results as compared with both DFT calculations p4j and simulations 



. In Ref. |39|] a fit to simulation results of the surface tension of a one-component hard- 



sphere fluid at a planar hard wall is presented, which gives quasi-exact results and closely 
resembles the SPT expression. 

In terms of the weighted bulk densities n , . . . ,n 3 , defined in Eqs. (p~5|) - (P~8|) , the SPT 
approximation for the surface tension of a hard-sphere mixture close to a planar hard wall 
241 can be written as 



n\ 1 ni 
PlsPT = ^- + ^~ n * ■ (26) 



i 

l-n 3 ' 87r(l-n 3 ) 2 ' 

This expression reduces to the one-component SPT approximation of the surface tension 
when it is evaluated for the one-component bulk weighted densities. In this latter case 
the surface tension can be expressed solely in terms of the packing fraction rj of this single 
component. 

Within the Rosenfeld functional the surface tension 7 of a binary hard-sphere mixture 
at a planar hard wall follows from the equilibrium density profiles Pi(z), i = s,b, as obtained 
in the previous subsection: 

pn + [3 P v 



A 



lim 

L— >oo 



(3pL + J dz U[ Ps (z),p b (z)] + J2i 

n I i=s,b 



(27) 



In contrast to the strong dependence of the results of the excess adsorption on the choice 
of functional, we find that the original Rosenfeld functional as well as its modifications 
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predict very similar results for the surface tension in the whole range of packing fractions 
studied here. Our results for the surface tension of a binary hard-sphere mixture with size 
ratio Rh : R s — 3 : 1, calculated within the original Rosenfeld functional, are shown in 



Fig. [H]. The deviation between these results and those obtained by the modifications of the 
Rosenfeld functional are at most 3% and deviate from the predictions of the SPT at most 
by 10%. 



IV. SUMMARY AND CONCLUSIONS 

Based on the Rosenfeld density functional we have analyzed structural and thermody- 
namic properties of binary hard-sphere mixtures near a hard wall with the following main 
results: 

1. Figures [L] and f| demonstrate that the structure of both density profiles Pi(z), i = s,b, 
of a binary hard-sphere mixture close to a planar hard wall as obtained by the original 
Rosenfeld functional is in excellent agreement with simulation results for size ratios 
Rb : R s = 5 : 3 and Rf, : R s — 3 : 1, respectively. The high level of agreement between 
our DFT results and the simulation data from Ref. fH| is confirmed quantitatively 
by small mean square deviations E s and Ej, defined in Eq. (|22|) . In terms of these 
quantities all versions of the Rosenfeld functional are practically of the same quality. 
Only at high packing fractions rather small differences between the results obtained 
by different versions of the Rosenfeld DFT become visible like those shown in Figs. |3| 
and [5|. 

2. The concentration profiles (see Fig. ^ and Fig. |J) calculated from the density profiles 
confirm the depletion picture: the small spheres are depleted from regions close to the 
wall while the big spheres are enriched. 

3. The numerical accuracy of our calculations is demonstrated in Tables | and [IT] by the 
high degree at which the sum rule Eq. (^), which relates the sum of the contact values 
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of both density profiles with the equation of state, is respected. The sum rule, however, 
makes no prediction for the individual contact values and we find in Tables [III] and [IV] 
that each version of the Rosenfeld functional takes a different route to satisfy the sum 
rule. 

4. Using the modified Rosenfeld functional corresponding to Eq.(Q) we have calculated 
the excess adsorption of the small spheres T s (rj s , r] b ) (see Fig. 0) and of the big spheres 
Tb(?7s, r)b) (see Fig. |8|) as functions of the packing fractions rj s and r\\, for a binary hard- 
sphere mixture with size ratio Rb : R s = 3 : 1. We find that these quantities depend 
very sensitively on the accuracy of the numerical calculations and, as can be seen in 
Fig. H they differ significantly from the excess adsorption calculated by the original 
Rosenfeld functional. 

5. The surface tension of a binary hard-sphere mixture with size ratio Rb : R s = 3 : 1 



close to a planar hard wall is shown in Fig. [It]. We find that all versions of the 
Rosenfeld functional give results which are in good agreement with the prediction of 
scaled particle theory [Eq. (pSD]. 

From these results we conclude that the class of Rosenfeld functionals yields quantita- 
tively reliable descriptions of interfacial structures in binary hard-sphere fluids. We expect 
that the same level of reliability also holds for multi-component hard-sphere fluids. 

The excess adsorptions r s and Tb of the small and big spheres emphasize the differences 
between the various versions of the Rosenfeld functional most. In order to decide whether 
the original Rosenfeld functional or whether its modifications predict these quantities more 
accurately, additional simulation data of the excess adsorption in a binary hard-sphere fluid 
are needed. 
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APPENDIX A: CALCULATION OF THE WEIGHTED DENSITIES 

Within the minimization procedure of the Rosenfeld functional the weighted densities 
n a and n a have to be calculated repeatedly. Therefore it is necessary to optimize these 
calculations with respect to both computational speed and numerical accuracy. 

The weight functions of the Rosenfeld functional are given by 

4 3 \v) = Q(\r\-R l ), (Al) 

4 2 \r) = 6(\r\-R z ), (A2) 

and 

«S 2) (r) = j^ ( 5(|r|- J R i ) (A3) 

with the Heaviside function G and the Dirac delta function 5. The remaining scalar weight 
functions are proportional to oof^: uo^p = /(47rRi) and uof^ = uf^ /(AuR^). The first 
vector weight function is collinear with lj^: oj^p—ujf ^ /(AnRi). 

In order to calculate the weighted densities integrals of the type 

1^ = J dV A (rK(r - r') (A4) 

have to be evaluated. For these convolution type integrals one can exploit the symmetry 
properties of the density profiles. For the present geometry the weighted densities can be 
written as 

= e r dz 'p^ + (as) 
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with s and b for small and big, respectively, and with reduced weight functions uo^ which 
are functions of z only: 

J?\z) = v(R\ - z\ (A6) 
uP(z) = 2nR i , (A7) 

and 

af t) (z) = 2nze g (A8) 

with the unit vector e 2 in z-direction. The relations between these and the remaining weight 
functions are the same as for the original weight functions. The integrals in Eq. (|A5|) are 
one-dimensional convolutions which can be calculated faster and more accurate in Fourier 
space than in real space. By introducing the Fourier transforms of the density profiles, 

Pi(k) = FT{ Pi {z)), (A9) 

and those of the weight functions, 

&<r\k)=jv(Q<r\z)), (aio) 

the weighted densities can be expressed as 



^rM^AlW)). (ah) 

\i=s,b / 



This route of calculation offers the important advantage that the numerical calculation of 
these convolutions can be speed up significantly by applying Fast-Fourier- Transform (FFT) 
methods. Moreover it turns out that calculations of convolutions in real space depend more 
sensitively on the grid size Az to be used for discretization than those in Fourier space. We 
expect that the reason for this is that the FFT algorithm interpolates between data points 
with trigonometrical functions. To overcome this problem in real space a sophisticated 
integration scheme would have be applied or a very small grid size would have to be chosen. 
Both remedies additionally slow down the numerical calculation in real space. 



The results presented in this appendix are applicable if the density profiles depend on 
the z coordinate only. However, similar results can be obtained if the density profiles have 
radial symmetry. 
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FIG. 1. The density profiles of the small spheres (a) and the big spheres (b) of a binary 
hard-sphere mixture with size ratio Rb : R s = 5 : 3 and packing fractions r/ s = 0.0607 and 
rjb = 0.3105 close to a planar hard wall (dashed line). a s = 2R S is the diameter of the small 
spheres. The solid lines correspond to results obtained by a free minimization of the original 
Rosenfeld functional and the symbols (□) denote simulation data from Ref. |i~0|| . The full dots 
indicate the values at contact. The positions of the extrema of p s (z) and pb(z) are very close to 
each other if they are measured from the corresponding position of contact, i.e., z = R s = ^a s 
and z = Rb = |cr s , respectively. The dotted lines correspond to the coarse grained density profiles 



defined in Eq. (23). They indicate that in spite of the high contact values the wall actually leads 
to a slight net depletion for both species in the first layer. 
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FIG. 2. The concentration profiles <& s (z) an d &b( z ) of the small spheres and big spheres, respec- 
tively, corresponding to the density profiles in Fig. [fl. For geometric reasons pi(z < Ri) = with 
R s = \o s and Rb = |crt,. Therefore $ s and §b are defined only for z > a s and &b(Rs < z < Rb) = 
and $ b (R s < z < R b ) = 1. 
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FIG. 3. In order to resolve the small differences in the density profiles of the small (a) and big (b) 
spheres between the DFT results obtained by the various versions of the Rosenfeld functional, parts 
of the density profiles shown in Fig. [j] are magnified here. The solid lines correspond to the original 
functional, whereas the dotted, dashed, and dashed-dotted lines correspond the interpolated version 
(int RF), and the antisymmetrized modification with q = 2 and q = 3, respectively. The parameters 
are the same as in Fig. ||. 
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FIG. 4. The density profiles of the small spheres (a) and of the big spheres (b) of a binary 
hard-sphere mixture with size ratio Rb : R s = 3 : 1 and packing fractions r/ s = 0.0047 and 
rjb = 0.3859 close to a planar hard wall. The solid lines correspond to results obtained by a free 
minimization of the original Rosenfeld functional. The simulation data (□) are taken from Ref. . 
As in Fig. [I] the positions of the extrema of p s (z) and pb(z) are very close to each other if they 
are measured from the corresponding position of contact, i.e., z = R s = ^a s and z = Rb = |cr s , 
respectively. The dotted lines correspond to the coarse grained density profiles defined in Eq. (|23|). 
Different from Fig. [I], here the net density of both species is clearly enhanced near the wall. 
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FIG. 5. In order to highlight the small differences of the DFT results obtained by different 
versions of the Rosenfeld functional, here small parts of the density profiles shown in Fig. || are 
magnified. 
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FIG. 6. The concentration profiles & s ( z ) and <&b( z ) of the small spheres and big spheres, re- 
spectively, corresponding to the density profiles of Fig. |j. R s = ^a s and Rb = |cr s so that due 
to geometric constraints &b(Rs < z < Rb) = and Q S {R S < z < Rb) = 1. The anti-correlated 
behavior of the two profiles is strongly locked in. 
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the packing fractions r\ s and rjb in the fluid phase |36|. The dots (•) denote turning points. 
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FIG. 8. The excess adsorption of the big spheres near a planar hard wall as function of the 
packing fractions rj s and rjb in the fluid phase []36|]. The dots (•) denote turning points. The square 
symbols (□) denote simulation results for the excess adsorption of a one-component hard-sphere 
fluid taken from Ref. |35j . 
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FIG. 9. Excess adsorptions T s and Tf, calculated with the original Rosenfeld functional com- 
pared with those obtained from its modification corresponding to Eq. (||). The packing fraction of 
the small spheres is r/ s = 0.15. 
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FIG. 10. The surface tension /?7(r? s , r?;,) of a binary hard-sphere mixture in the fluid phase |3C] 
with size ratio Rf, : R s = 3 : 1 as function of the packing fractions rj s and calculated from 
the original Rosenfeld functional (lines). For r] s =0 the hard-sphere fluid freezes at 776=0.494. The 



square symbols (□) denote simulation results from Ref. [38] for the one-component hard-sphere 
fluid. 
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TABLES 

TABLE I. Contact values at a planar hard wall of a binary hard-sphere mixture with size 
ratio Rb : R s = 5 : 3 for various packing fractions rj s and The comparison of the sum of 
the contact values with the Percus-Yevick compressibility equation of state, which underlies the 
Rosenfeld functional, tests the accuracy of our numerical procedures. For small total packing 
fractions rj = rj s + 775 this equation of state is in good agreement with the more accurate equation 
of state Ppmcsl established by Mansoori et al. [j3l|. However, for larger values of rj, there are 
deviations. 
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TABLE II. Results for the same quantities as in Table | for a binary hard-sphere mixture with 
size ratio Rj, : R s = 3 : 1. 



V 


Vs 


Vb 


l^iPi{Ri + o) 


Q C 

PPpy 


PPMCSL 


0.1209 


0.0188 


0.1021 
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0.1288 


/~l AAT7 
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r\ r\r\ o -i o 
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0.0283 
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0.2230 


0.0026 
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0.00643 
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0.00637 


0.2471 


0.0199 


0.2272 


0.01590 


0.01590 


0.01578 


0.2513 


0.0089 


0.2424 


0.01104 


0.0H04 


0.01094 


0.3021 


0.0016 


0.3005 


0.01180 


0.01180 


0.01159 


0.3257 


0.0136 


0.3121 


0.02166 


0.02166 


0.02124 


0.3775 


0.0099 


0.3676 


0.02837 


0.02837 


0.02747 


0.3906 


0.0047 


0.3859 


0.02718 


0.02719 
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TABLE III. Individual contact values of the density profiles for a binary mixture with size 
ratio Rb : R s = 5 : 3 as obtained by the original Rosenfeld functional (afp° r9 (z = Ri + 0)), the 
antisymmetrized modification with q = 3 (afpf~ 3 (z = R4 + 0)), and the interpolating modification 
{afpf lt (z = Ri + 0)) with i = s,b. Oi = 2Ri is the diameter of species i. While all versions of 
the functional respect the sum rule in Eq. (pQ), the individual contact values differ slightly for the 
different versions of the DFT. 77 = r] s + iji, is the total packing fraction. 



V 


* 3 sPs r9 (Rs + ) 


a 3 sP r 3 (R s +) 


alpf t (R s +) 






4p? t (Rb+, 


0.1246 


0.07051 


0.07058 


0.07059 


0.34182 


0.34150 


0.34144 


0.1252 


0.27332 


0.27340 


0.27343 


0.14522 


0.14487 
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0.38657 


0.47860 


0.47619 
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0.13385 


0.13432 


0.13444 


0.82502 


0.82284 
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0.2247 


0.05002 


0.05036 


0.05043 


1.13782 


1.13628 


1.13594 


0.2278 


0.90612 


0.90702 


0.90727 


0.36210 


0.35787 


0.35677 


0.2670 


1.15640 


1.15871 


1.15943 


0.70205 


0.69099 


0.68829 


0.2749 


0.27715 


0.27944 


0.28026 


1.67448 


1.66398 


1.65994 


0.3477 


0.17094 


0.17520 


0.17709 


3.43107 


3.41124 


3.40174 


0.3712 


0.47247 


0.48279 


0.48746 


4.07483 


4.02673 


4.00460 


0.3906 


0.03099 


0.03160 


0.03181 


2.39204 


2.38927 


2.38826 


0.3968 


0.10016 


0.10471 


0.10703 


5.33002 


5.30870 


5.29677 
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TABLE IV. Results for the same quantities as in Table | for a binary mixture with size ratio 
Rb '■ Rs = 3 : 1. 



rj a 3 s p° s r9 (Rs+) <T 3 s p q s =3 (Rs+) <T 3 s Pl nt (Rs+) *ip° b r9 (Rb+) crlpl= 3 (R b +) a 3 pT\R b +) 



0.1209 


0.04862 


0.04864 


0.04864 


0.36808 


0.36769 


0.36761 


0.1288 


0.00932 


0.00933 


0.00932 


0.42401 


0.42387 


0.42389 


0.1456 


0.07891 


0.07895 


0.07895 


0.50360 


0.50248 


0.50231 


0.2230 


0.00839 


0.00842 


0.00841 


1.16319 


1.16241 


1.16248 


0.2471 


0.07074 


0.07101 


0.07104 


1.52497 


1.51744 


1.51686 


0.2513 


0.03121 


0.03136 


0.03136 


1.54220 


1.53814 


1.53814 


0.3021 


0.00626 


0.00633 


0.00633 


2.37903 


2.37714 


2.37714 


0.3257 


0.05996 


0.06072 


0.06078 


3.05954 


3.03862 


3.03744 


0.3775 


0.05168 


0.05289 


0.05299 


4.73342 


4.69975 


4.69719 


0.3906 


0.02540 


0.02612 


0.02617 


5.18598 


5.16579 


5.16450 
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